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We argue that colossal magnetoresistance is a critical phe- 
nomenon and propose a mechanism to describe it. The mech- 
anism relies on the halfmetallic behavior of the materials 
showing colossal magnetoresistance, and yields a correlated 
percolation model that, we argue, captures all qualitative fea- 
tures of colossal magnetoresistance, above as well as below 
the Curie temperature. The model only serves for revealing 
the underlying mechanism of colossal magnetoresistance, and 
does not aim to reproduce precise, numerical results. 
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I. INTRODUCTION 

In recent years, much experimental effort has been de- 
voted to materials displaying Colossal Magnetoresistance 
(CMR), which is a strong dependence of the resistance on 
the magnetic field as well as on temperature. The studied 
materials are rare earth manganese perovskites |lp|)| , but 
recently the effect has been shown |l(J to occur also in 
T^M^Oy, which has a pyrochlore structure. These ma- 
terials show a ferromagnctic-to-paramagnetic transition 
at a certain temperature, the Curie temperature. The 
resistance of the materials shows a peak at or near this 
temperature, and falls off quite rapidly with higher and 
lower temperatures. Switching on a magnetic field also 
lowers the resistance. It is this behavior of the resistance 
as a function of temperature and magnetic field that is 
called CMR. 

The manganese perovskites are of a mixed valence 
type; if only Mn 3+ or Mn 4+ are present the material 
is insulating and antiferromagnetic, and no CMR is ob- 
served. It is known for a long time that the double ex- 
change interaction between pairs of Mn 3+ and Mn 4+ is 
responsible for the ferromagnetic and metallic properties 
of the perovskites p| . In this picture, a dependence of 
the conductance on the spin direction of the charge car- 
riers is already present. More recently band structure 
calculations ]l2^ | using the local spin-density approxima- 
tion indicated an effective halfmetallic behavior, and in 
similar calculations adopting the generalized gradient ap- 
proximation the halfmetallic character fully emerged [jl3| . 
The double exchange mechanism thus is responsible for 
the halfmetallic behavior of these materials, be it that 
the valence electrons are less localized than in the origi- 
nal model of references 

It is observed in all measurements of CMR that the 
peak in the resistance occurs close to or at the Curie 



temperature T c . At this temperature the spontaneous 
magnetization of the materials vanishes; it is a critical 
point. The occurrence of the peak in the resistance at or 
close to this temperature is a strong indication that CMR 
itself is a critical phenomenon, that is, that the behavior 
of the resistance is intimately connected with the critical 
behavior of the domains of magnetization. If this is true, 
typical band structure calculations require, due to the 
infinite correlation length, infeasible large system sizes 
to reproduce the effect. 

The goal of this paper is to propose a mechanism for 
CMR. Relying on the halfmetallic character of the mate- 
rials displaying CMR, we will introduce a model that, we 
believe, captures the basic features of the phenomenon, 
but does not aim to reproduce the correct numerical val- 
ues of relevant temperatures, resistances and so on. In- 
deed, our proposed model is much too simple as com- 
pared to the no doubt very complicated processes that 
govern CMR. On the other hand, revealing the basic 
mechanism that is responsible for the typical features of 
CMR is the first step towards a satisfying understanding 
of the phenomenon. 

Moreover, in our proposed mechanism, CMR is a crit- 
ical phenomenon. The concept of universality in critical 
phenomena then assures that the universal features of 
the phenomenon, such as qualitative behavior, scaling 
functions, and critical exponents, will be reproduced cor- 
rectly, as they do not depend on the precise definition of 
the model, but only on general features like symmetries 
and dimensionality. Other quantities, like the precise lo- 
cation and height of the peak in the resistance, do depend 
on the precise definition of the model, and will, regarded 
the extreme simplicity of our model, be reproduced in- 
correctly. The universal properties, however, can be put 
to strict experimental tests to decide upon the validity of 
the proposed mechanism. 

The mechanism we have in mind is that of correlated 
percolation, and the typical model that brings in the re- 
sistance is the correlated resistor network. The model is 
introduced in section |l| In section III we will shortly ex- 
plain the concepts of percolation and resistor networks, 
and discuss the phase diagram of our model. In this sec- 
tion, we will argue that the qualitative features of CMR 
that are exhibited in experiments are captured by the 
model. In section |y] we present Monte Carlo calcula- 
tions on the resistor network that make our results more 
explicit. In section ^ we will discuss several complica- 
tions that arise when the model is extended to become a 
more sophisticated explanation of CMR. We end with a 
conclusion. 
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II. A MODEL FOR COLOSSAL 
MAGNETORESISTANCE 

Materials displaying CMR turn out to be halfmctallic 
p|-^3|, that is, the band structure of the electrons de- 
pends on the relative orientation of the spin with respect 
to the magnetization of the material. Electrons with spin 
parallel to the local magnetization are conducting, elec- 
trons with spin antiparallel are insulating. The conduc- 
tance of the material thus depends on the structure of 
the domains with different directions of magnetization. 
As suggested in Ref. [||, this explains the observation 
that the resistance increases with temperature, as more 
domain walls occur with higher temperatures, thereby 
hampering the percolation of charge carriers. It is less 
clear that this effect also explains the decrease in resis- 
tance above T c , but we will see that it does, so there 
is no need to invoke another mechanism (e.g., magnetic 
polarons) to explain the behavior above T c . 

The simplest conceivable model that captures such a 
mechanism is a lattice model, where on every lattice point 
a vector m is defined, which represents locally the di- 
rection of magnetization in the material. Due to the 
extended wave function of the electrons, such a lattice 
point does not have to be identified with one atom; the 
lattice can be defined on a more coarse grained level as 
well. The direction of magnetization m must in some 
way be governed by an energetic interaction between the 
different domains, in such a way that for temperatures 
above the Curie temperature T c the average magnetiza- 
tion vanishes, whereas below the Curie temperature there 
is a finite magnetization. The allowed orientations of m 
in each domain must follow the magnetic anisotropy of 
the material, which is in the case of the CMR materials 
a perovskite-like symmetry. 

It is not our goal to introduce a precise, quantitative 
model that describes the perovskites. We want to have 
a simple model that only reveals the basic mechanism 
of CMR. Let us therefore, for simplicity, assume that m 
can only point in two different directions, up and down. 
We assume, then, that m obeys an Ising symmetry, so 
let us replace m by the more familiar notation of a scalar 
spin S € {+1,-1}. The interaction between the different 
domains is also defined in the most simple manner, that 
is, there exists a 'nearest neighbor' coupling K between 
the spins. In addition, we allow for an external magnetic 
field h. The Hamiltonian governing the statistics of the 
model is in this case the familiar Ising Hamiltonian 

H=-Kj2SiS j -hJ2Sj, (1) 

(ij) 3 

where (ij) denotes a summation over nearest neighbor 
pairs only. 

The resistance is brought into the model following the 
half-metallic character of the CMR materials. Resis- 
tances are defined independently for charge carriers with 



spin up and down. Charge carriers with spin up have a 
low resistance in domains characterized by S = +1 and 
a high resistance in domains with S = — 1. For these 
spin-up carriers we define a unit resistance on bonds be- 
tween two sites having S = +1 and an infinite resistance 

on (H — ) and ( ) bonds. The infinite resistance on a 

( ) bond reflects the insulating character and is a sim- 
plification, as these domains will have a large but finite 
resistance. The infinite resistance on (-1 — ) bonds means 
that our model does not allow for spin-flip processes, and 
this is no doubt a simplification as well. 

For the charge carriers with spin down a similar assign- 
ment of resistances is made. These assignments yield an 
expectation value of the overall resistance of the mate- 
rial independently for up and down carriers. We denote 
the inverse overall resistance, i.e., the conductance, by 
(<t + ) and (er_) respectively. The total overall conduc- 
tance then simply is 

(a> = (a + ) + (<7_). 

The angular brackets denote Ising expectation values. 

The assignment of resistances for the spin up carriers is 
shown in figure [I], where sites with S = +1 are depicted 
black, and sites with S = — I are white. In the figure 
resistances are defined between nearest as well as next- 
nearest neighbor sites; the reason for this is explained 
below. 




FIG. 1. The assignment of resistors between neighbor and 
next-neighbor sites for the charge carriers with spin up. The 
small circles denote the lattice sites; black sites have a magne- 
tization S = +1, white sites have S = — 1. If both sites have 
S = +1 a unit resistance is put on the bond between them, 
otherwise the bond receives an infinite resistance. A similar 
assignment is made for the charge carriers with spin down. 

The conductance of the model is calculated in the usual 
way: let T be a configuration of magnetizations Sj on the 
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lattice sites. The conductance of this configuration can 
be calculated from the assignment of resistances, and is 
denoted by er(r). The expectation value of the overall 
conductance then is obtained by a weighted summation 
over all configurations in the usual way, 



|^>(r) exp(-/3i/(r)), 



where (3=1/ kT, Z is the partition function of the Ising 
model and H(T) is the Hamiltonian of equation (|l|). We 
will now see what the predictions of this model are for 
the expectation value of the overall conductance of the 
lattice. 



III. THE PERCOLATION PHASE DIAGRAM 

The overall conductance of the lattice is defined as fol- 
lows: consider a square lattice of L x L sites with resis- 
tances on bonds defined as in figure [l]. Keep the lower 
row of the lattice at a fixed potential V = and the up- 
per row at V — 1 . The conductance of the lattice is then 
equal to the expectation value of the current. 

Conductance over the lattice is possible (e.g., non-zero) 
only if there exists a path from border to border exclu- 
sively over resistances R — 1, that is, if at least one of the 
two spin directions is percolating |[14 15|. The qualitative 
features of the model thus can be obtained by consider- 
ing the percolation phase diagram of the Ising model. 
In percolation, one defines clusters by putting bonds be- 
tween neighbor spins S = +1 or S = — 1. Let us consider 
percolation for sites S = +1; percolation for S = — 1 
is of course just the same but with the Ising magnetic 
field h replaced by —h. In our case, a bond is placed 
between each pair of nearest and next-nearest neighbor 
sites having S = +1. These bonds make up clusters, and 
the percolation order parameter P is the density of sites 
in the percolating cluster. If there is percolation then P 
is finite, if there is no percolation, P = Q. 

Our case is called correlated percolation as the dis- 
tribution of percolating (S — +1) and non-percolating 
(S = —1) sites is a correlated one. When this distri- 
bution is that of the Ising model, the phase diagram is 
known, and this is one of the reasons of choosing the Ising 
model as an example of our model: many exact results, 
also for percolation, are known for this model. Its per- 
colation phase diagram is known ||l6f when percolating 
bonds are placed between each nearest neighbor pair of 
sites having S = +1. In another paper [p.7[ , we derive 
the phase diagram when next-nearest neighbor percolat- 
ing bonds are defined as well. It is this phase diagram 
we need for our explanation; it is shown in figure @. 



percolating 



non-percolating 



FIG. 2. The percolation phase diagram for the Ising model 
for percolation of the up spins. Percolation is possible between 
nearest and next-nearest neighbor up spins. T is temperature 
and h is the magnetic field. The thick, solid line is a critical 
percolation line. The order parameter (the fraction of sites 
present in the largest cluster) is finite in the shaded area and 
goes continuously to zero when approaching the critical line. 
The critical line merges smoothly with the T-axis at T c , which 
is a tricritical point for percolation. The dashed line is a 
first order transition for percolation, meaning that the order 
parameter jumps from a finite value to zero on this line. 

In experiments, CMR occurs in thin films j|-|5| or 
even in bulk samples [|l],[3]|J. Our model thus should 
be pseudo three-dimensional or fully three-dimensional. 
This is the reason that in our two-dimensional model 
we defined resistances on next-nearest neighbor bonds as 
well. These 'crossing' bonds make the bond-graph non- 
planar, and this is the simplest way to mimic pseudo 
three-dimensional behavior. By defining these crossing 
bonds we allow for more possible paths between perco- 
lating sites, and the same is happening in a system con- 
sisting of several layers. 

Inclusion of these crossing bonds is crucial, because 
without the crossing bonds there is no percolation |Ttj| ] 
above the Curie temperature. At first it seems therefore 
that the percolation mechanism cannot explain the de- 
creasing resistance above T c . Inclusion of the crossing 
bonds, however, makes the region above T c percolating, 
and this is necessary for our explanation. 

The effect of even more possible paths or of bonds with 
an even larger percolation range is, however, quite small. 
Even in a fully three-dimensional model, where there are 
many more possible paths for a cluster to follow, the 
percolation point lies only a few percent below the Curie 
temperature We expect, therefore, our conclusions 
also to hold in the case of a different lattice model or a 
different dimensionality. The percolation point may lie 
somewhat below T c but will be close to it. The phase 
diagram in that case is shown in figure pi 
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the T-axis where no percolation occurs for both spin di- 
rections is the Curie point T c , so there is a peak in the 
resistance at T c . (In fact, in our model this peak is in- 
finitely high due to the choice R = oo on the insulating 
bonds.) All other points on the T-axis are percolating 
and hence show a better conductance. 

2) The asymmetric shape of the peak in the resistance 



FIG. 3. The phase diagram for a lattice model where the 
percolation point T p and the Curie point T c do not coincide. 
In this case the percolation temperature is always lower than 
the Curie temperature. The thick solid line is a critical per- 
colation line in the universality class of random percolation 
and the same is expected for the endpoint of the line. The 
dashed line remains a first order transition line. 

We choose to follow figure ^| for the explanation of 
our mechanism and the comparison with experiments. 
The qualitative comparison, however, holds for the other 
phase diagram of figure [| as well. 

For T — y oo or K — ► the probability distribution 
becomes a random one and the type of percolation is 
random percolation. The thick line merging smoothly 
with the T-axis is in the universality class of random 
percolation: on this line the order parameter vanishes 
algebraically with the same exponent as that for ran- 
dom percolation. The Curie point T c is the critical point 
of the Ising model, where there is a ferromagnetic-to- 
paramagnetic transition. For percolation, this point is a 
tricritical percolation point. The dashed line is a first or- 
der transition for the magnetization, and for percolation 
as well; the order parameter P jumps from a finite value 
to zero over this line. 

In principle, the percolation phase diagram tells noth- 
ing about the value of the conductance, only whether it 
is zero or non-zero. Nevertheless we may well assume 
that the closer one moves to the percolation threshold, 
the lower the conductance will be. From this, we can 
extract the qualitative behavior of the overall resistance 
from the phase diagram and draw a comparison with the 
experimental properties of CMR. To check the predic- 
tions of the phase diagram, we performed Monte Carlo 
calculations on the Ising resistor network. The results 
of these calculations confirm the predictions of the phase 
diagram, and are discussed in the next section. 

Our conclusion is that the qualitative features of CMR 
are all present in the phase diagram of correlated perco- 
lation. Consider the following features: 

1) The peak in the resistance at the Curie temperature 
T c . According to the phase diagram, the only point on 



as a function of T. This shape is clearly present |QM| 10 
in experiments. In the phase diagram, moving away from 
T c is moving away from the percolation point and hence 
decreasing the resistance. Due to the presence of the 
critical percolation line that merges smoothly jl6| with 
the T-axis, we expect this decrease to be slower above T c , 
resulting in an asymmetric peak in the resistance. This 
expectation is confirmed using a Monte Carlo calculation, 
described in the next section. 

3) The peak in the resistance as a function of mag- 
netic field. In every experiment, the resistance peaks at 
zero field, a feature that is clearly present in the phase 
diagram. Moving away from the zero field axis means 
stimulating the percolation of one of the spin directions 
and hence decreasing the resistance. Note that our model 
does not give a correct numerical value of the magne- 
toresistance ratio AR/R, where AR is the difference in 
resistance with and without magnetic field. What does 
follow from the model, however, is the qualitative result 
that this ratio is largest at T c . 

4) A particular observation in experiments is the shape 
of the peak in the resistance as a function of the mag- 
netic field. As can be seen, e.g., from Fig. 1 in Ref. 0, 
Fig. 5 in Ref. ||], and Fig. 5 in Ref. ||, plots of the re- 
sistance versus the field display a cusp at temperatures 
below T c but are smooth at higher temperatures. This 
remarkable feature follows directly from the percolation 
phase diagram: below T c there is a first order transition 
line. Moving away from this line in both directions of the 
field immediately decreases the resistance, meaning that 
a plot of the resistance versus the field shows a cusp, the 
peak being strongest at T c . Above T c , on the other hand, 
there is no transition line at zero field, and hence the be- 
havior of the resistance as a function of the magnetic field 
must be smooth. 

5) It is observed §§110 in CMR that the P eak in 
the resistance shifts to higher temperatures when a mag- 
netic field is switched on. This directly follows from the 
presence of the percolation line in the phase diagram at 
temperatures above T c . Switching on the magnetic field 
above T c one 'feels' the vicinity of this line, whereas below 
T c the nearest percolation point is T c itself, which is fur- 
ther away. Hence the resistance decreases more strongly 

function of the field below T c . 
This comparison shows that the basic features of CMR 
that are displayed in all experiments, are directly ex- 
plained by our simple percolation mechanism. One single 
and simple mechanism accounts for its features, above as 
well as below the Curie temperature. 
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IV. RESISTANCE CALCULATIONS 

To make some predictions of our model, as explained in 
the previous section, more explicit, we performed Monte 
Carlo calculations on the Ising correlated resistor net- 
work. The model we used is that described in section || 
with the Hamiltonian of equation (|l|) and the assignment 
of resistors of figure [l]. 

We performed Monte Carlo simulations of the expec- 
tation value of the conductance for several lattice sizes. 
We used the Wolff-algorithm [jl9| for the Monte Carlo 
part, and the multigrid method of Edwards, Goodman 
and Sokal [^0|, based on the standard code AMG1r4 [ fjl| , 
to calculate the conductance of a spin configuration. The 
latter code is slightly changed, in order to cope with re- 
sistors on next-nearest neighbor bonds. 

The interest in the first place is the behavior of the 
resistance for different temperatures. We performed cal- 
culations on different system sizes; the result is shown 
in figure H The peak in the resistance at T c becomes, 
for an infinite system, infinitely high, a result already 
following from the percolation phase diagram. The cal- 
culations confirm also the asymmetric shape of the peak 
in the resistance. 
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FIG. 4. Resistance as a function of temperature for differ- 
ent system sizes at zero field. The system sizes are indicated 
in the plot, the resistance is scaled such that the minimal 
resistance (all bonds having R = 1) is unity. The lines are 
guides to the eye. Resistors are defined as in figure [lj For 
larger systems, the peak shifts to T c , becomes infinitely high 
but remains asymmetric. For T <T C there is a rapid conver- 
gence to the minimal (unit) resistance. 

Secondly, the interest is in the critical exponents of the 
correlated resistor network. The interested reader is re- 
ferred to reference (l7) for a more detailed account. In 
this reference, we describe our calculations of the conduc- 
tance exponent t of the correlated resistor network. This 



exponent governs the algebraic decay of the conductance 
upon approaching a percolation threshold. Two different 
exponents play a role in the phase diagram of figure |^. 
The first is the exponent of the critical percolation line. 
It is defined as 



a(h) 



\h - h c 



where <x(/i) is the conductance at fixed T > T c , and h c is 
the critical value of the magnetic field. The critical per- 
colation line is in the universality class of random perco- 
lation, so the value of the exponent t on this line must be 
that of the random resistor network. The random resistor 
network is a long standing unsolved problem in statisti- 
cal physics, but good numerical results of the exponent 
t are known. The best estimate [£2| in two dimensions 
known to us is t = 1.299 ± 0.002. 

The second exponent involved is that governing the 
vanishing conductance at the Ising critical point. For 
percolation, this point is in a different universality class, 
and therefore the conductance exponent also differs from 
the random one. We calculated this exponent from the fi- 
nite size behavior of the conductance, as explained in ref- 
erence [0 . Because the Ising critical point is a tricritical 
point for percolation, there are two non-equivalent direc- 
tions of approaching the critical point with two different 
exponents involved. For the temperature direction, the 
exponent t is defined as 



a(T) ~ \T-T C \* atfc = 0, 



(2) 



yielding a value t = 0.2000 ± 0.0007. For the other, field- 
like, direction, the exponent is defined as 

a(h) ~ \h\ l for T = T c . 

In this case, the exponent t has the value t = 0.1067 ± 
0.0004. 

As will be explained in the next section, the numerical 
values of these exponents arc of no direct interest, as the 
Ising model is not the correct lattice model for describ- 
ing the CMR-materials. What is of interest, however, is 
the difference between the exponent of the critical line 
(random resistor network) and those of the Ising criti- 
cal point (correlated resistor network). These exponents 
differ considerably. Even when there turns out to be no 
tricritical point in the correct phase diagram of the CMR- 
materials, as in figure ||, the model will, in a sense, be 
close to the tricritical point, which means that there will 
be an observable crossover from the random exponent to 
the value of the tricritical exponents. 



V. DISCUSSION 

We chose the Ising model to serve as an example for our 
mechanism for CMR, because of its simplicity and for the 
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fact that many exact results are known for this model. 
There are, however, several features of our model that 
render it unrealistic for a more sophisticated explanation 
of CMR. The first is of course that the Ising model is not 
the appropriate lattice model, as the Ising variables Sj 
do not have the symmetry of the CMR-materials. The 
latter are mostly perovskites, so the appropriate lattice 
model should have a magnetization vector mj at each site 
j that follows the symmetry of the perovskites. 

A further shortcoming of our model is the infinite re- 
sistance for spin-up charge carriers on (H — ) and ( ) 

bonds. In principle, we should allow for a finite (but 
large) resistance on these bonds. The effect of this can be 
compared with the ferromagnetic-to-paramagnetic phase 
transition in the Ising model: if we switch on a small 
magnetic field, the true phase transition disappears, but 
the magnetization, specific heat and so on still show signs 
of the vicinity of the critical point. The same will happen 
when introducing these large but finite resistances. The 
actual value of this resistance displays the possibility of 
a spin-flip mechanism of the charge carriers. The impor- 
tance of this effect is, we believe, quite small, such that 
the true phase transition is smeared out but the behavior 
of the overall resistance still shows signs of the vicinity 
of the percolation transition. 

The finite height of the peak in the measurements can 
simply be incorporated by a large but finite 'background' 
resistance i?^ parallel to the net resistance of the corre- 
lated resistor network, representing the resistance of the 
'non-conducting' bonds. The total resistance as a func- 
tion of temperature T, following from equation (^|), then 
becomes 



1 + all — In 



(3) 



where T p is the temperature where the resistance peaks, 
so T p is equal to or slightly lower than the Curie temper- 
ature T c . The amplitude a is non-universal and has in 
general different values above and below T p . 

Experiments so far have not been set up to measure 
critical exponents. For the exponent t of equation (j^), 
a rough analysis of the published data seems to point to 
a value between 1.5 and 2.5. It would be worthwhile to 
have more accurate experimental data on the exponents. 

On the theoretical side, the interest is in the critical 
exponents of the correlated resistor network for the spin 
model corresponding to the magnetic anisotropy of the 
perovskites. Calculating these exponents is considerably 
more elaborate than in the Ising case, as even the perco- 
lation phase diagram has, to our knowledge, never been 
studied for models other than Ising. 



VI. CONCLUSIONS 

We presented a correlated percolation mechanism for 
CMR and argued that its observed effects can be ex- 
plained by considering the materials displaying CMR as 
a correlated resistor network. The introduced model only 
sheds light on the underlying mechanism, but is certainly 
not the appropriate one to reproduce accurate numeri- 
cal results. For this, the model has to be changed into 
a spin model displaying the right magnetic anisotropy 
of the perovskites and has to be made pseudo or fully 
three-dimensional. Nevertheless, we argue, the qualita- 
tive results of the model will remain untouched. 
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